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We investigate black-body induced energy shifts for low-lying levels of atomic systems, with a 
special emphasis on transitions used in current and planned high-precision experiments on atomic 
hydrogen and ionized helium. Fine-structure and Lamb-shift induced black-body shifts are found 
to increase with the square of the nuclear charge number, whereas black-body shifts due to virtual 
transitions decrease with increasing nuclear charge as the fourth power of the nuclear charge. We 
also investigate the decay width acquired by the ground state of atomic hydrogen, due to interaction 
with black-body photons. The corresponding width is due to an instability against excitation to 
higher excited atomic levels, and due to black-body induced ionization. These effects limit the 
lifetime of even the most fundamental, a priori absolutely stable, "asymptotic" state of atomic 
theory, namely the ground state of atomic hydrogen. 

PACS numbers: 12.20.Ds, 31.15.-p, 42.50.Hz 



I. INTRODUCTION 

Energy shifts of atomic levels due to interactions with 
black-body radiation are primarily important for Ryd- 
berg states, where large dipole polarizabilities enhance 
the magnitude of the effect [3, 0| , but at the current 
rate of advance of high-precision spectroscopy, its impor- 
tance for the correct realization of frequency standards 
has also been stressed in the literature 0, [H, IE 0j- ^ n 
particular, the work of Farley and Wing [2J was mainly 
focused on Rydberg states, which were accessible to high- 
precision spectroscopy at the time. In contrast, we here 
focus on the ground state and on the 2S state of hydro- 
gen and hydrogenlike ions. The transition frequencies 
in these systems have only recently come within reach of 
high-precision laser spectroscopy, due to the development 
of frequency combs. In all cases studied by Farley and 
Wing Q, the Lamb shift and fine structure shifts could 
be neglected, but this does not hold for the transitions 
studied here. 

Ionization by black-body radiation also is an important 
effect for excited atomic states. The ionization process 
leads to a finite width of the states (resonances), with 
the energy acquiring a (small) imaginary part. The sign 
of the imaginary part is negative, and we can write the 
black-body induced energy shift A£bb as 

A£ bb = AE hh - i^J± (i) 

where A£"bb is the real part of the energy shift, and 
AFbb/ft is the black-body induced width. 

Several peculiarities characterize the black-body in- 
duced radiative shift. First of all, we recall that two vir- 
tual processes contribute: one where the atom absorbs 
a black-body photon and then returns to the reference 
state by emission into the same mode of the electromag- 
netic field from which a photon had been absorbed, and 
another one where the sequence of absorption and emis- 



sion processes is reversed. This is very much analogous 
to the so-called ac Stark shift (see, e.g., [8|, [9() that an 
atom feels in a laser field, but with the difference that 
the black-body radiation is isotropic, and that the polar- 
ization vectors are equally distributed among all possi- 
ble directions in space. For the ground state being the 
reference state, the "width" (imaginary part of the en- 
ergy shift) is exclusively generated by the absorption-first 
channel, where one of the propagator denominators be- 
comes resonant and the virtual state takes the role of 
the final state of the ionization or excitation process [l(| • 
Note that recently [ll[ , a finite width has been predicted 
for the ground state of hydrogen at very high tempera- 
tures on the basis of field-theoretical considerations [see 
Eq. (23) of Ref. [Tl|]. 

One observation, potentially of fundamental interest, 
substantiated here by a concrete calculation, is that even 
the IS state of atomic hydrogen, acquires a finite width 
when the hydrogen atom is exposed to black-body radi- 
ation, although the ground state of the most fundamen- 
tal atomic system is otherwise assumed to be the perfect 
"asymptotic state" used in all S'-matrix type calculations 
regarding the quantum electrodynamic (QED) shifts of 
energy levels [H, EH G3- However, before we come to 
an analysis of this effect, we briefly revisit the evaluation 
of the real part of the black-body induced energy shifts, 
for the 1S-2S transition in hydrogen [IE [IE] and ion- 
ized helium (see Sec. [II]). We then analyze the imaginary 
parts of the black-body induced energy shifts in Sec. IIHI 
Conclusions are drawn in Sec. IIVI 



II. REAL PARTS OF BLACK-BODY 
RADIATION INDUCED ENERGY SHIFTS 

For clarity, we keep all factors of h, c and eo in all 
calculations. The Boltzmann constant is denoted ks- We 
recall that the energy distribution per frequency interval 
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dv of black-body radiation is 

, , , 8nhv 3 
p{v) dv — - — exp 



hv 
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which is connected to the time-averaged square of the 
electric field e 2 (v) at frequency v (in units of the square 
of the electric field strength per frequency) as 

e 2 (v)dv = — p(v)dv. (3) 

We write the Cartesian coordinates as x l and rewrite the 
dipole polarizability of an nS state to take the isotropic 
character of the polarization vectors of black-body radi- 
ation into account [cf. Eq. (20) of Ref. [9(]. Then, the 
black-body energy shift AEbh{nS) is given by the for- 
mulas 
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H - E(nS) ± hv 



iS 



(4a) 



TABLE I: Black-body energy shift of the ground state 
of atomic hydrogen and ionized atomic helium, at a 
temperature of T = 4 K, T = 77 K and T = 300 K. The 
shift A_Ebb(lS l ) is obtained by numerical integration ac- 
cording to Eq. (I4b|l and divided by the Planck constant 
in order to be expressed in frequency units. 



Nucl. charge Temperature 
number 


h AE bh 


(IS) 


Z = 1 


4K 


-1.22 x 10 


" 9 Hz 


Z = 1 


77 K 


-1.68 x 10 


" 4 Hz 


Z = 1 


300 K 


-3.88 x 10 


" 2 Hz 


Z = 2 


4K 


-7.65 x 10 


- 12 Hz 


Z = 2 


77 K 


-1.05 x 10 


" 5 Hz 


Z = 2 


300 K 


-2.42 x 10 


3 Hz 



oo 

AE hh (nS) =- — (P.V.) / dve 2 {v)P„(nS) 



(4b) 



where E(nS) is the energy of the reference nS state. The 
latter integral actually has to be taken as a principal- 
value (P.V.) integral because the propagator denomina- 
tors in (|4"aj) can become singular when a black-body pho- 
ton hits an atomic resonance (see the discussion below in 

secnni. 

For the IS state, the properties of the dipole polar- 
izability P V (1S) are well understood [13, and P V (1S) 
can be evaluated using an entirely nonrelativistic ap- 
proximation [IH, and H in (|4a| may be replaced by the 
Schrodinger Hamiltonian. To see why that is the case, 
let us consider the typical photon frequencies at which 
the black-body photon energy peaks. Ignoring the "—1" 
in the denominator of (J3j , it is easy to derive the result 



3k B T 



(5) 



for the frequency at which p{y) assumes its maxi- 
mum value. We consider a temperature range from 
4K...300K. In this temperature range, the frequency 
i/ max varies from 2.5 x 10 11 Hz to 1.9 x 10 13 Hz, which is 
well below the Rydberg constant expressed in frequency 
units, i?oo c = 3.298 x 10 15 Hz, which defines the fre- 
quency range for the virtual excitations to P levels that 
are decisive for the evaluation of the dipole polarizability. 
For the ground state acting as the reference state, we may 
therefore even ignore the terms ±hv in the propagator de- 
nominators in (|4a|) and replace the dynamic by the static 
polarizability, which can be evaluated easily and is used 
here in the nonrelativistic approximation (cf. Ref. [l||) 



Po(lS) = I 
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(6) 



The dipole polarizability defined in Eq. (|4a[) parametri- 
cally scales as Z~ 4 , because the energy differences in the 
propagator denominators scale with Z 2 , whereas the two 
dipole transition matrix elements in the numerator each 
acquire a factor 1/Z (the atom becomes "smaller" by a 
factor Z^ 1 as the nuclear charge number increases). We 
finally obtain the approximation 



AE hh (lS) 
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The corresponding data in Table U are in excellent agree- 
ment with this approximation. 



A subtle poin t concerns the dipole polarizability of ex- 
cited states [19(: While for the 15 state, it is entirely suffi- 
cient to approximate H by the Schrodinger Hamiltonian, 
this is not the case for the 2S state and for higher excited 
nS states with n > 2. Indeed, the 2P 1 / 2 state is displaced 
from 2S only by the Lamb shift, and the 2P 3 / 2 state is 
displaced in energy only by the fine-structure, and yet, 
the dipole transition matrix elements between the 2S and 
the 2P states are manifestly nonvanishing. The dipole 
transition matrix elements which enter the numerators 
for the dipole polarizability, can still be taken in a non- 
relativistic approximation, but the denominators must be 
adjusted for the Lamb shift and the fine-structure effects. 
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In summary, we have for an nS level with n > 2, 
P v {nS) 



1 



1/2 



EE E • 
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\(nS 1 / 2 (m=±) \x z \ nP 1/2 (m = fi))\ 
E(nP 1/2 ) - E(nS) ± /«/ 
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E 
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EE 
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E(nP 3/2 ) - E(nS) ± hv 



\lltt E 



\{nS(m = Q) \x l \ = 
£(n'P) - E(nS) ± hv 



(8) 



Here, the \nS\/ 2 ), the |nP!/ 2 ) and the \nP 3 / 2 ) states 
can be approximated by Schrodinger-Pauli wave func- 
tions, i.e. the radial part is taken in the nonrelativis- 
tic approximation, but the angular part is given by a 
spinor (two-component) function xjl(^) defined as in [2fJ; 
k = 2 (l—j) (j + 1/2) is the Dirac angular quantum num- 
ber, with I being the orbital angular momentum quantum 
number and j the total electron angular momentum, and 
[i E { — (j + 1/2), . . . , j + 1/2} is the magnetic projec- 
tion (half- integer). The states nP and 2S, by contrast, 
are plain nonrelativistic Schrodinger eigenstates with an 
angular part of the form Y/ m (f), where I is the orbital 
angular momentum, and m e {— I, . . . , 1} is the magnetic 
projection (integer). 

In the following, we concentrate on the case n = 2, 
i.e. on the black-body shift of the 2S state. The key 
to finding an analytic approximation to A_Ebb(2S*) is the 
following. We observe that for most cases of interest 
studied here, the condition 



E fs -C h v max < h Roo c , 



(9) 



is fulfilled, where c « 3.289 x 10 15 Hz is the Rydberg 
constant expressed in frequency units. This hierarchy 
means that we can make different approximations for the 
first two as opposed to the third term on the right-hand 
side of ([5]). Namely, the relevant frequency range of the 
black-body radiation is large as compared to the fine- 
structure and the Lamb shift, but small in comparison to 
the frequencies corresponding to transition with a change 
in the principal quantum number. For the last term, 
of ([H]), we can thus make the static approximation and 
obtain 



wv= i-EEE E 



= 1 /j=-1 n'=3 



= 120 



1 \(2S(m = 0)\x i \n'P(m = 
3 E{n'P) - E(2S) ± hv 

h \ 2 1 



m e c J (Za) A m e c 2 



(10) 



where we reemphasize that the sum over n! starts from 
ml = 3. The corresponding energy shift is 



A^ b (2S) 



16tt 3 4 T 4 
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(11) 



For the black-body shift corresponding to the first two 
terms on the right-hand side of ([5]) , we first recall the fine- 
structure interval and the Lamb shift in leading order, 



E fs = E(2P 3/2 ) - E(2S 1/2 



(Za) 4 m c c 2 
Z2 



E L = E(2S 1/2 ) - E(2P 1/2 ) = — (Za) 4 m e c 2 
x 1 1 In [(Za)- 2 ] + ^ (lnfc (2P) - lnk (2S)) + 



(12a) 
(12b) 
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Here, we make the opposite approximation and expand 
the terms for large v, 
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Using ([3]), we can integrate integrate over the black- 
body spectrum and obtain 



AK b (25) = 



5m c cr 



a 



' B 



Am r c 2 



[\ln[(Za)- 2 } 
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+ - (lnfco(2P)-lnfco(25)) + -| T 2 Z 2 . (14) 

The approximation then is 

AP bb (25) w A£' bb (25) + AP b ' b (2^) . (15) 

From Table [TT] we see that the approximation holds for 
all cases studied except for the shift of the 2S state for 
ionized helium (Z = 2) at 4K. The inequality f9| is 
not fulfilled in this case, because z/ max = 2.5 x 10 11 Hz 
at 4K whereas E {s /h = 1.8 x 10 n Hz for Z = 2 (see 
also Fig. [T]). For Z = 1, at the same temperature, the 
inequality Ef s <C h v max is better fulfilled, and this affords 
an explanation for the fact that the approximation works 
well in the latter case. 

One would naively assume that the black-body shifts 
should decrease with the nuclear charge, because of the 
Z~ 4 scaling of the polarizability. However, that is not 
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FIG. 1: The black-body spectrum given by (|2)| is displayed as a function of the frequency v and of the temperature T. In order 
to visualize the spectrum over a wide range of temperatures, we plot p 1 ^ 10 , and we use arbitrary units (a.u.), normalized so 
that p assumes a value of unity at its maximum (as a function of v) for T = 400 K. The maximum of the blackbody spectrum 
as a function of the frequency at given temperature is denoted as v max [see also Eq. J5}] . The relation of v max to the relevant 
atomic frequencies, as immediately discernible from the plot, illustrates why the approximation and (|14[) is not applicable 
to ionized helium at T = 4K. The frequency vis-?,p is the frequency for excitation of the 2S state to the 3P state for ionized 
helium. 



the case. The formulas (ITT|) and (| 14[) indicate that there 
are two competing effects for the shift of the 2S state, one 
which is related to virtual transitions with a change in 
the principal quantum number [see Eq. ((lip ] . and another 
one which is related to fine-structure and Lamb shift 
transitions [see Eq. ([14]) ]. The former scales with T 4 Z~ 4 , 
as expected, but the latter scales as T 2 Z 2 , and increases 
with the nuclear charge number (somewhat counterintu- 
itively). We would thus like to refer to this as an anoma- 
lous scaling. This particular behavior, in connection with 
planned experiments (2lT |. has been a major motivation 
for carrying out the calculation report here. 

Incidentally, we note that the black-body shifts at 
T = 300 K are by several orders of magnitude larger than 
other non-standard, non-resonant effects ("accuracy lim- 
its") for two-photon spectroscopy recently discussed in 
Refs. [iH H3l- Note that after leaving the cooled nozzle 
in the current hydrogen atomic beam experiment [H|,[lB|, 
the slow hydrogen atoms enter a high vacuum which is 
kept at room temperature. 



III. IMAGINARY PART OF BLACK-BODY 
INDUCED ENERGY SHIFT 

In order to evaluate the black-body induced decay 
rates, we have to be more precise than in Eqs. (I4a|) and 
(I4b[) . Namely, we have to collect all the poles from the 
bound-state spectrum as well as integrate over the transi- 



tions to the continuous spectrum, and introduce infinites- 
imal imaginary displacement into to propagator denom- 
inators. The relevant formulas are 



P„(nS) 



i=l 



1 



H - ie - E(nS) - hv 
1 

H -ie- E(nS) + hv" 



nS 

nS)) , (16a) 



AS hh (nS) = - lim %- f Ave 2 {v)P v {nS) 

e^0+ 2 J 


= AE hh (nS) - i . 



(16b) 



The first term on the right-hand side of (|16a[) describes a 
process with an absorption of a photon by the atom from 
a black-body mode, which lowers the energy of the black- 
body photon field in the virtual state (term —hv), and 
a subsequent emission into a black-body mode. For this 
term, the imaginary part is generated for virtual atomic 
states, as contained in the spectral decomposition of H, 
which have an energy higher than that of the reference 
state E(nS). This is the only relevant process to gener- 
ate a black-body induced width if the ground state acts 
as a reference state, and it is easy to see that the imagi- 
nary part of the energy shift thus generated has the cor- 
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TABLE II: Black-body energy shift of the 2S excited state of atomic hydrogen and ionized atomic helium, evaluated 
according to Eq. (JSJ), expressed in terms of frequencies via division by the Planck constant h. Note that the entries 
for Ai5 bb (2»S') + Aii' b ' b (2,S') are given here only as an indication of the quality of the approximation given in Eqs. (JTTJ 
and (|14[1 and therefore in brackets. The entries for AE bb (2S) are obtained by numerical integration of (I4b[l and 
thus relevant for comparison to experiment. The explanation for the discrepancy of the approximation and the 
numerical integration for 4 K and Z — 2 is given in the text (see also Fig. [1} . 

Nucl. charge number Temperature h' 1 AE bb (2S) h' 1 [AE bb (2S) + AE bb (2S)] 



z = 


1 


4K 


7.79 x 


10" 


7 Hz 


(8.13 x 


10" 


- 7 Hz) 


z = 


1 


77 K 


-1.44 x 


10" 


3 Hz 


(-1.46 x 


10" 


3 Hz) 


z = 


1 


300 K 


-9.89 x 


10" 


'Hz 


(-9.87 x 


10" 


'Hz) 


z = 


2 


4K 


3.40 x 


10" 


6 Hz 


(3.30 x 


10" 


■ 5 Hz) 


z = 


2 


77 K 


1.18 x 


10" 


2 Hz 


(1.19 x 


10" 


2 Hz) 


z = 


2 


300 K 


1.21 x 


10" 


'Hz 


(1.21 x 


10" 


"'Hz) 



TABLE III: Black-body induced width AF bh (lS)/h of 
the ground state and of the 25* excited state, evaluated 
according to Eq. (|16b|) for atomic hydrogen. The results 
are given in Hz (cycles per second). 



State Temperature h 1 AF bb (nS) 



1,5 


300 K 


1.13 x 


10" 


163 


1,5 


3000 K 


2.15 x 


10" 


9 a 


IS 


30000 K 


8.00 x 


10 6 


Hz 


2,5' 


300 K 


3.08 x 


10" 


3 H: 


2,5 


3000 K 


7.49 x 


10 3 


Hz 


2,5' 


30000 K 


2.34 x 


10 7 


Hz 



rect sign, i.e. it contributes a negative imaginary part to 
A£ hh (nS). 

The second term on the right-hand side of (|16ap de- 
scribes a process with (first) emission into a black-body 
mode, then absorption from the black-body field. In this 
case, the imaginary part is generated for virtual atomic 
states of an energy lower than that of the reference state 
E(nS). This process can be relevant for the generation 
of a black-body induced decay when the 2S states acts as 
a reference state, and the virtual state is a 2P\/2 state. 
It is again easy to sec that the imaginary part of the en- 
ergy shift thus generated has the correct sign, i.e. it also 
contributes a negative imaginary part to AS^^nS) . 

We can then easily obtain the black-body induced de- 
cay rates of the IS* and 25* states. Numerical results are 
given in Table IIII1 where we restrict ourselves to atomic 
hydrogen. There is a very strong increase of the rates 
with temperature, and we note a comparatively large in- 
duced decay rate for the 2S state, due to the proximity 
of the 2P\/2 and 2P 3 / 2 levels. 



IV. CONCLUSIONS 

In this paper, we have set up a formalism for the 
treatment of black-body induced energy shifts and corre- 
sponding decay rates (widths) for low-lying atomic states 
in atomic systems with a low nuclear charge. We have 
evaluated our expressions for the 15* and 25* states of 
atomic hydrogen and hydrogenlike (ionized) helium. The 
formula ([8]) leads to a consistent treatment of the virtual 
nPj states (j — 3/2 and j = 1/2) which are displaced 
from the reference nSi/ 2 state only by a fine-structure 
splitting and by the Lamb shift, respectively (see Sec. HI)) , 
and can be easily generalized to other cases of interest. 
We note that all black-body shifts given in Tables III and ITU 
are well below 1 Hz in absolute magnitude and do not 
have to be taken into account at current and projected 
levels of accuracy for high-precision laser spectroscopy, 
although it is perhaps useful to remark that they are 
larger than other non-standard effects ("accuracy lim- 
its") for two-photon spectroscopy recently discussed in 
the literature (Refs. HHH). 

The imaginary part of the black-body induced energy 
shift ("black-body ac Stark shift") has been discussed 
and evaluated in Sec. lIIIl The corresponding decay rates 
given in Table UTTl can be quite substantial, especially at 
elevated temperatures. At T — 3000 K, the 2S state 
of hydrogen acquires a black-body width of about 8 kHz, 
and at T — 30000 K, the ground-state of atomic hydrogen 
is roughly 8 MHz wide, purely due to the black-body in- 
duced interactions. At the same temperature, the Boltz- 
mann factor exp[— [E^p — E\s) / {ksT)] for the excitation 
of the ground state into the 2P state is only 2%. 

Although the black-body energy shifts have been stud- 
ied quite intensively, the numerical results and approxi- 
mations have not yet appeared in the literature to the 
best of our knowledge. We have worked in SI units 
throughout this article in order to enhance the clarity 
of the details of the derivation. 
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